Collective modes of spin-orbit coupled Fermi gases in the repulsive regime 
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We investigate the collective modes in the spin-orbit coupled Fermi gas with repulsive s-wave 
interaction. The interplay between spin-orbit coupling and atom-atom interactions plays the crucial 
role in the collective behaviors of Fermi gas. In contrast with ordinary Fermi liquid, spin-orbit 
coupled Fermi gas has strongly correlated spin and density excitations. Within the scheme of random 
phase approximation, we classify collective modes based on the symmetry group and determine their 
properties via the poles of corresponding correlation functions. Besides, the particle-hole continuum 
is obtained, where the imaginary part of these correlation functions become non- vanishing. We make 
comparisons with collective excitations in the ordinary Fermi liquid without spin-orbit coupling and 
in a helical liquid, i.e., surface states of a three dimensional topological insulator. We also propose an 
experimental protocol for detecting these collective modes and discuss corresponding experimental 
signatures in the ultracold Fermi gases experiment. 
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I. INTRODUCTION 

A great deal of attentions have been focused on the 
spin-orbit coupling (SOC) because of its fundamental in- 
terests in condensed matter systems [lHl] and important 
applications in spintronic device Q- In recent years, a 
wide range of atomic physics and quantum optics tech- 
nology provides unprecedented manipulation of a variety 
of intriguing quantum phenomena, therefore it seems to 
provide an ideal platform to study the effects of SOC in 
ultracold atomic systems. The experimental studies on 
this topic have made great breakthrough. Based on the 
Berry phase effect Q and its non-Abelian generalization 
Q , Spielman's group in NIST has successfully generated 
a synthetic external Abelian or non-Abelian gauge po- 
tential coupled to neutral atoms [Iol - tl3j . Recently, the 
SOC Fermi gas has been first engineered in weakly re- 
pulsive K [14[ or 6 Li |l5[ atomic gases. Realization of 
SOC in quantum gases will open a whole new avenue in 
cold atom physics. 

Motivated by these recent experimental progresses of 
ultracold Fermi gases, we consider the two dimensional 
(2D) Fermi gas with SOC in the repulsive regime. The 
repulsive atom-atom interaction can be engineered in the 
upper branch of energy spectrum, where there are uncon- 
densed Fermi gas in the absence of molecule formation 
(lfjj . The repulsive Fermi gas is stable when it is far 
away from the resonant regime, which has been success- 
fully reached in experiment [17]. In our previous work, 
we have studied low energy single particle excitations and 
calculated the Fermi liquid parameters such as the quasi- 
particle lifetime, renormalization factor and the effective 
mass in the repulsive regime fl8j . 

In this paper, we investigate the collective modes of two 
dimensional Fermi gas with SOC in the repulsive regime. 



The research of the low energy collective modes in degen- 
erated quantum gases yields a wealth of insights into the 
properties of ultracold atomic systems. A lot of previous 
studies have devoted to the collective behaviors in vari- 
ous symmetry broken phases of ultracold atomic systems, 
which include the degenerate gas at the Bose-Einstein 
condensation (BEC) to Bardeen-Cooper-Schrieffer (BCS) 
crossover [l9|-l22j , Fermi gas in the unitarity limit [231 - 4261 ) , 
imbalanced Fermi gas [27j, and BEC in the presence of 
SOC [H,[2§|. By contrast, we focus on the normal state 
regime of the SOC Fermi gas in this work. The interplay 
between SOC and s-wave atom-atom interactions plays 
the crucial role in our investigation. Compared with or- 
dinary Fermi gas, SOC Fermi gas has strongly correlated 
spin and density excitations. Within the scheme of ran- 
dom phase approximation (RPA), we classify all the col- 
lective modes and determine their properties via the poles 
of corresponding RPA correlation functions. 

In contrast with previous works in solid state sys- 
tems, the consideration of s-wave interaction in ultra- 
cold atomic systems instead of the Coulomb interaction 
leads to qualitatively different collective behaviors. The 
reasons for this are twofold: First, the force range of 
the s-wave interaction is short and the interaction vertex 
is independent to the momentum transfer. Second, the 
s-wave interaction is spin-dependent in ultracold Fermi 
gas, which can be decomposed into the density and spin 
channels [3(|. The collective modes of this system are 
grouped into two categories: (i) one branch of gapless 
mode, namely zero sound, which is an oscillation of den- 
sity coupled with the transverse spin oscillation; (ii) three 
branches of gapped modes. We expect our microscopic 
calculations of the collective modes would have imme- 
diate applications to the SOC Fermi gas in the upper 
branch of the energy spectrum. 

The paper is organized as follows. The model building 
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of SOC Fermi gas with repulsive s-wave interaction is de- 
scribed in Sec. II, where all the microscopic parameters 
and the helicity eigenstates are explained. Furthermore, 
the collective modes are classified based on the symmetric 
property. In Sec. Ill, we develop a general approach to 
calculate a series of RPA correlations functions. In Sec. 
IV, we investigate the solutions of matrix forms of RPA 
equation analytically and numerically. In Sec. V, we pro- 
vide a comparison and discussion of the collective modes 
in various fermionic systems. Finally, we propose an ex- 
periment protocol to detect collective modes in Fermi gas 
with SOC on the upper branch of energy spectrum and 
estimate corresponding experimental signatures in Sec. 
VI. 
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II. THE SOC FERMI GAS WITH REPULSIVE 
S-WAVE SCATTERING 

We consider the collective modes of a 2D spin-1/2 ul- 
tracold repulsive Fermi gas with Rashba SOC, described 
by the model Hamiltonian, 

U =XH,A^ c k,/3 + 2 5 ^ 4+ q ,T c P-q,i c P^ c M-> (1) 

k,a,/3 k,p,q 

where a and /3 are the spin indexes. The hrst part is the 
non-interacting Hamiltonian, and h is the single particle 
Hamiltonian with SOC 

k 2 

h = h A(z X k) • a - fj,, (2) 

2m 

where /i = k F /2m is the chemical potential and kp is the 
Fermi momentum in the absence of SOC, A represents 
the strength of SOC. Here and bellow, h is taken as 1. 
The eigenstates of single particle Hamiltonian h can be 
obtained as follows 

|k,±l) = ^( ± J, (k) ), (3) 

where </>(k) = arctan(fc a /fc x ) is the azimuthal angle of 
momentum k and the helicity ±1 represents that the 
in-plane spin polarization is right-handed or left-handed 
with respect to the momentum. Hence the eigenstates 
in Eq. ([3]) are naturally called helicity bases. The 
dispersion relations for two helicity bases are £k,± = 
(k 2 ± 2fcfl|k| — k F )/2m, where Ur = m\ corresponds to 
the recoil momentum in experiments [lil . [l5j . The Fermi 
momentum and the recoil momentum provide two scales 
for the 2D Fermi gas with SOC, and the dimensionless 
ratio of them 7 = k^/kp denotes the significance of SOC. 
The Fermi surfaces are given by £ks = 0, which gives rise 
to two circles in the momentum space with two Fermi 
momenta as k s — nkp — ska, where k = ^Jl + j 2 . We 
plot the energy spectrum and the Fermi surfaces with 
the associated spin textures in Fig. [T] for two different 
chemical potentials. We note that the outer and inner 



FIG. 1. (Color online) (a) and (c) plot the energy spectrum 
in presence of spin-orbit coupling with different fillings. The 
thick black horizonal line denotes the level of chemical poten- 
tial, (b) and (d) show the Fermi surfaces and the associated 
spin textures corresponding to (a) and (c) respectively. 



Fermi surfaces shown in Fig. Q] (b) have two opposite 
helicity, and therefore one has Berry phase ir, another 
— 7r. While the outer and inner Fermi surfaces shown in 
Fig. [T](d) have the same helicity, and therefore the Berry 
phase 7r and 7r respectively. The system experiences Lif- 
shitz transitions when the chemical potential passes the 
Dirac point of the spectrum at k = [3l|, [32j- In the 
present work, we only consider the case shown in Fig. [1] 
(a), where the two Fermi surfaces correspond to different 
helicity bands. 

The second term in Eq. (JT|) represents the s-wave 
interaction in ultracold atomic gases. The low energy 
interaction among ultracold atoms is universally deter- 
mined by the scattering length a s [33-35]. For quasi-2D 
system, which can be realized through a strong confine- 
ment in the z direction perpendicular to the 2D plane, 
the effective s-wave scattering strength is determined by 
2g = 4TrNa s /3^/2~7Tm( z , where N is the total atom num- 
ber, m is the mass for atom, and ( z = \J I /mui z is the 
confinement length of the atomic gases along i direction 
with uj z as the trap frequency of the confinement poten- 
tial. The significance of s-wave interaction in ultracold 
atoms could be characterized by the dimensionless ra- 
tio of the average interaction energy and kinetic energy 
Sint/zkin, which equals rag = 2irNa s /3V2nm( z . With 
the technique of Feshbach resonance, the strength of the 
repulsi ve interaction can be tuned within a wide range 

The dynamical response of SOC Fermi gas can be de- 
scribed via the density and spin susceptibility 

X ^ (q, iu m ) = E E T dre^(T T ci a a^c k+clP (T) 

k,p a/3jS Jo 

xc P+q7 ^c p5 (Q)), (4) 



where = (<7°, <r). We could study the collective modes 
through the poles of the dynamical response function to 
external density and spin perturbations. The SOC is 
isotropic in the 2D x-y plane. The density and spin sus- 
ceptibility is invariant under the simultaneous rota- 
tion of the momentum q and the spin s = \cr around the 
z axis, namely the rotation group: SO(2)x>, where the 
subscript T> represents the combined operation for the 
momentum and spin space. 

Furthermore, considering the Rashba-type SOC and 
the rotation invariance of Fermi surface around the z 
axis, the density excitations are coupled with the trans- 
verse components of the spin excitations s T = (z x q) • s 
intrinsically, while both of them are decoupled from the 
longitudinal com pon ent s L = qs and perpendicular com- 
ponent s z = z-s [43]. Therefore, it is convenient to study 
the collective modes under the helical spin bases, which 
is defined as {n, s T , s L , s z } with n as the density compo- 
nent. Under the helical spin bases, the 4x4 susceptibility 
matrix decomposes into two 2x2 matrices. One is in the 
subspace of n-s T , and the other is in the subspace of 
s L -s z . Based on this decomposition of the susceptibil- 
ity, the collective modes are grouped into two categories: 
one with the density and transverse spin excitations s T , 
and the others with the perpendicular and longitudinal 
spin excitations s z ,s L . Without loss of generality, we 
can assume q = ge x in the following. Thus we have 



(a) 



(b) 



8 T = 8V,8 L = 8* 



(5) 



In the following sections, we will work in this representa- 
tion to calculate the density and spin susceptibility and 
investigate the collective modes of SOC Fermi gas. 



III. DENSITY AND SPIN CORRELATION 
FUNCTIONS 



A. Feynman Rules and Susceptibilities 



The particle line shown in Fig. [5] (a) represents the 
Green's function with SOC in the Matsubara formalism 



G° (k,ife„) = 



,=±i lkr > 



(0) 



where ik n = (2n + l)7rfcsT is the fermionic Matsubara 
frequency and P s is the projection operator to the helicity 
eigenstates: P s = [1 + s(z x k) ■ <x]/2. It is helpful to 
apply the Fierz identity 2e^e Q ^ = (a^f 1 (ct m ) q5 to the 
s-wave interaction vertex, where ct m = (cr° , <x) and <r M = 
(<7°, — <r). The interacting part can be rewritten as 

V = 2T2T E E VaS,p-r4 (k + q) c\ (p - q) 



C/3 (P) C Q (k) , 



P, 1P„ 
* 



k+q,y p-q,5 

q 




(c) 



p+q. ip„ + iq„, 



(d) 



© = o + OO 

+ ^^^A/W^^^/VW^^^ + 



FIG. 2. (Color online) (a) The free Green's func- 
tion with SOC: — G°(p, ip„). (b) The interaction vertex 
— g (cr M )' 97 (u M ) c *' 5 . (c) One bubble diagram: the bare sus- 
ceptibility x A "'( < l! *9m). (d) The sum of the bubble diagrams 
gives the RPA susceptibility XiipAili iQm)- 



where V aS ,p 7 = g (c M )^ 7 (<7n) aS - The Feynman rule for 
the interaction vertex is shown in Fig. [5] (b). In fact, 
the interaction vertex considered here includes two cases: 
(k f,p |) -> (k' t,p' I) and (k t,P I) -> (k' i,p' t). 
It's important to notice that the infinite summation with 
this interaction vertex actually include both bubble and 
ladder diagrams [43| . 

The collective modes and particle-hole pairs are two 
fundamental types of excitations of the SOC Fermi gas. 
The dispersions of collective modes are determined by the 
poles of the density and spin susceptibility. With these 
Feynman rules defined above, the RPA susceptibility can 
be evaluated as (see Fig. d](d)) 

X % PA (q, iq m ) = X W (q, (— V~ ^) P "> (») 

1 + mx (q, nm) 

where 



V = diag{l, -1,-1,-1}, 



(9) 



(7) 



and x is the bare susceptibility for noninteracting Hamil- 
tonian. 

In Matsubara formalism, the bare susceptibility 
X Miy (q, iq m ) is given by the one loop diagram as shown in 
Fig. [1(c) 

X^ (q, iq m ) = -k B T J2 Tr[G° (k + q/2, ip n ) 
k,ifc n 

x G° (k- q/2, ik n -iq m )a„], (10) 

where iq m = 2m7rfcsT and ik n = (2n + l)7rfcsT are 
the bosonic and fermionic Matsubara frequencies, re- 
spectively. The Green's function G° (k, ik n ) includes a 
momentum-dependent projection operator P s , and the 
trace in Eq. (jTTJJ) gives rise to an overlap factor as 

F£ - Tr [P. (k + q/2) o»P r (k - q/2) a v \ , (11) 

where /x, v = 0,1,2,3 represent the density and spin 
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FIG. 3. (Color online) Particle-hole continuum of SOC Fermi 
gas for 7 = 0.5. The red region surrounded by the thick 
black lines represents the interband particle-hole continuum. 
The region surrounded by the dashed yellow lines represents 
continuum of intraband particle-hole excitations with helicity 
+1, the blue region filled with vertical lines, with helicity — 1. 
The points a, b, c and d correspond to the momenta 2ka, 2fe+i, 
2nkF and 2fc_i respectively, where the static susceptibility 
function exhibits singular behaviors [44|. 



components in x,y,z directions, and s,r are helicity in- 
dexes. After summing over the fermionic Matsubara fre- 
quency ik n , and performing the analytical continuation 
iq-rn — > cj + i0 + , we have 



/ (fk-q/2,a) - / (Ck+q/2,r 



£k-q/2,s - £k+q/2,r + W + i0+ 

. (12) 

At zero temperature, the Fermi occupation function is 
/(£) = 9(-£). The numerator in the Eq. (TljJ) is non- 
zero only if 



q| = 2fc_i caused by the intraband virtual transitions 
44j. These specific momenta are shown in Fig. G2 which 
appear as edges of the intra- and inter-band particle-hole 
excitation continuum. 

We want to make some general observations of the 
overlap factor Fj£ u and bare susceptibility x as follow- 
ing: (i) Ft. (q) = F sr (q). (ii) Fj? (q) = Ftf (-q). (hi) 
X (q, w) = % (— q, — w). In the long wave length limit 
|q| — > 0, the overlap factor F^ v can be expanded rel- 
ative to q up to 0(q 2 ) approximately. The intraband 
contribution is 



( 



\ 



1 

-s cos 9 
ssin# 



-scost/ 
cos 2 9 

sin 29 
2 

s|q| sin 29 



s sin f 

sin 29 
2 

sin 2 9 

■ s|q|sin 2 
1 2\kT 



■ |q| sm 9 

*-2|kT 

■ s|q| sin 26 

1 — ijk] — 

- s|q| sin' 6 
! 2|k| 

o 



/ 



(15) 

where 9 is the azimuthal angle of k = ( cos 9, sin 9) . The 
interband contribution is 



F^ v _ = 



( 



\ 





s|q| sin 2 9 

2|k| 
s|q| sin 29 

, 4|k| 

■ |q|sin_9 
1 2|k| 



2|k| 

sin 2 9 

sin 29 
2 

is sin 6 



»|q|« 



i 2rt 



4|k| 

sin 29 



is cost 



■ q sm ( 

•*~2|kT 
-is sin£ 

-is cos ( 

1 



(16) 

In the following calculations, we will find that the sus- 
ceptibility decomposes into two 2x2 matrices, which 
coincides with the argument in Sec. II. B based on the 
symmetry property of this system. Besides, the momen- 
tum q dependent terms in the overlap factor have negli- 
gible contributions to the low-g (long wave length limit) 
properties of the collective modes, such as the sound ve- 
locity of the gapless mode and the energy gaps of the 
gapped modes. 



Sk-q/2,s > 0> £k+q/2,r < 0, 



(13) 



Intraband contributions (r = s) 



Ck-q/2 )S < 0, 6c+q/2,r > 0, (14) 

which represent the conditions for particle-hole excita- 
tions. The £k-q/2,s — £k+q/2,r m the denominator is 
the corresponding particle-hole exciton energy, which in- 
cludes the contribution from intraband and interband 
and forms a continuum as shown in Fig. [3] Accordingly, 
the contributions to the susceptibility come from two as- 
pects: the intraband r — s and the interband r = — s. 
When the frequency and momentum of the susceptibility 
w,q fall in the particle-hole continuum (see Fig. [3]), the 
integration in Eq. (|12[) will develop a non-zero imaginary 
part, which corresponds to the damping of the collective 
excitations. The details about this integration will be 
further discussed in the following sections and Appendix 
A. The static density-density susceptibility x°°( t l) ex- 
hibits singular behaviors at |q| = 2fc# and |q| = 2nkp. 
In addition, there are weak anomalies at |q| = 2/c+i and 



In the long wave length and low frequency limit (q <C 
kp, uj <C /i), the expansions of the particle-hole exciton 
energy and occupation functions in Eq. (Ti"2"|) to the lead- 
ing order of |q| are 

c t k + sk R ~ 

£k-a s -£k+a, s ^ k-q, (17 

2 2 m 

and 

/(£k-f ,s) ~ /(&+§,*) - (6k,») (4-f ,s - &+§,«)■ 

(18) 

Both have corrections of 0(q 2 /kp). The intraband sus- 
ceptibility reads as 

„„ / \ m v-^ k s f ,„_.„, cos 9 

yT q, iq m ) = -T^~J2t~ / d0F ss— 5> 

4tt z k kp J y + —cos9 

(19) 
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where y = muj / nkp\q\ and 9 is the angle between k and 
q. Making use of the overlap factor given in Eq. (TT5j). 
the elements of x (q, w ) can be expressed in terms of the 
following integrals 



'(q,w) = 



/ h lh 

lh h o 

h o 

V o ooo 



(20) 



where the definitions and results of these integrations 
Iii I21 13, a n d I5 are given in Appendix A. 



C. Interband contributions (r = — s) 

Similar to the intraband contribution, the energy of the 
interband particle-hole excitation is expanded relative to 
q to the leading order 



m V 



2sfcfl — k ■ q 



(21) 



with corrections of 0(q 2 /k 2 F ). We notice that the en- 
ergy of interband particle-hole excitation obtains an ad- 
ditional term ±2kku/m in the presence of SOC. This 
term dominates the energy of particle-hole excitation at 
the region q <C &r, which is dramatically important for 
giving rise to the low-g properties of the collective exci- 
tations, including the sound speed of the gapless modes 
and the energy gaps of the gapped modes. The difference 
of occupation is expanded as 



/(£k-§, s ) - /(Ck+§,- 8 ) 
+[5(k-k a ) + 5(k-k- 



-sQ (k R 
,kq 



nkp\) 



(22) 



The difference of occupation contains two contributions, 
which give rise to interband susceptibility: ^ mtcr = 



^,mter _|_ ^-jn tcr , xhe second part is negligible for prop- 
erties of the collective modes in the low-g regime where 
|q| < k R . 

The first part xl nter is an integration within the mo- 
mentum shell nkp — k R < |k| < nkF + &r, which can be 
evaluated straightforwardly for weak SOC with 7 <C 1 



, inter 

Xi 



7r|q| 



E 



V 





2k R 






7|q|^4 











SKI2 —inli 

inli skIq 



\ 



(23) 



The definitions and results of these integrals I\, I2, I a 
are given in Appendix A, and the arguments of them 
are y + 2skn/\q\. The additional term 2sk R /\q\ comes 
from the split of the energy band. For strong SOC with 
7 ~ 0(1), the evaluation of ^J ntor are simplified in the 
long wave length regime with q <C k R 



inter 
Xi 



m x - 

IT 



/ 

F s /2 

F s /2 

Vo F s 



•O(^), (24) 



where F s = (27 - 



• z 



J±f±^)/8 7 with z s = W /4s 7 e F . It 
can be shown that Eq. (|24|) reduces to Eq. ((23)) to the 
first order of 7 at the low-g regime. 

The second part xii tcr (Qj w ) is an integration per- 
formed in two separated Fermi surfaces 



Xlf er (q,w) 



m x - ± 

7T ' Sir 



+ 



8tt 



) cost 



cos + z0+ 



where y s = ^ 



y s — cos 6 + i0+ 



(25) 



2s T'^> = " 2s q ■ Resembling 
the results of the intraband contribution ^ lntra (q, w), we 
can also express the results in terms of some well defined 
integrals 



, .inter 

Xn 








$[h(v.)-h(v.)] 
h(y s )} h{y s ) + h{y s ) 








h(y s ) + h(y s ) 
is[h (ys) - h (y s )} 



\ 


-is[I 2 (y s ) - h iVs)] 
h(y s ) + h(y s ) ) 



(26) 



All of these integrations and corresponding results are 
given in Appendix A. 

The summation of the intraband and interband con- 
tribution gives the total susceptibility 

X = X intra + Xi nt ° r + Xlf (27) 

For the low-g behaviors of the collective modes, we find 
that the contribution of Xn ter is less significant. This 



could be understood from that the second term of the 
right-hand side in Eq. (|22[) is much smaller than the first 
one in the low-g regime |q| <C k R . Therefore the energy 
gaps of the gapped modes and the sound speed of the 
gapless mode are mainly determined by ^ lntra and x™ ter . 
Furthermore, we will find that the gapped modes are 
mostly determined by the interband contribution Xi nter 
in the following section. 
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a =2.25a„ 




which are given by 



FIG. 4. (Color online) (a) The speed of zero sound as a 
function of s-wave scattering length, (b) The speed of zero 
sound as a function of the strength of SOC. The parame- 
ters used here are: the number of 40 K atoms is about 10 4 , 
k,R — 2ir/\ with A = 773nm, 7 = 0.5, trapping frequencies 
(u}±,u} z ) — 2ir x (10,400)Hz, and a s = 2.25ao, where ao is 
the Bohr Radius. The corresponding dimensionless interac- 
tion strength mg is about 2.5, which is less than the critical 
value 7r. 



IV. COLLECTIVE MODES 

In this section, we study the collective behaviors in the 
long wave length regime with |q| <C ka via the poles of 
the density and spin RPA correlation functions [46| . The 
susceptibility can be expanded with respect to the small 
dimensionless parameter |q|/fc_R. 

We consider the weak SOC at first, for which we have 
y s ~ y s ~ y + 2sfcfl/|q|. In the long wave length limit 
|q| <C kn, there is \y s \, \y s \ ^> 1. In Appendix A, we give 
the asymptotic behaviors of the integrals in Eq. (l20l) . 
(|23]l and p6|) for \y s \, \y s \ ^> 1. Given these asymptotic 
behaviors, the susceptibility in the subspace of n-s T can 
be expanded to O ((|q|/fc^) 2 ). To study the low-g be- 
haviors of the collective modes, it is useful to note the 
following behaviors of x^ v ■ If the dimensionless ratio 
y = mu I 'nkp |q| is fixed when q — > 0, we obtain 



XnT 



Xnn 



^ yXnn 
- + V 2 Xnn 



|q| ->■ 0,y fixed, (28) 



where \nn = —mli{y)/ir is the density-density element 
of the susceptibility matrix. The formula of \nT in the 
long wave length limit with y fixed holds for strong SOC. 
If we fix the w when |q| — > 0, the dimensionless parameter 
y tends to infinity, which leads to 



XnT 








|q| -> 0,w fixed, (29) 



where the intraband has no contribution in this limit. 
From Eq. (|28|) and ([29]) . we conclude that the solutions 
contain two different categories: (i) the gapless modes, 
which are coupled oscillations of the density and trans- 
verse spin excitations; (ii) the gapped modes, which are 
the oscillations of the transverse spin excitations. 

The collective modes in this subspace can be deter- 
mined by the poles of RPA susceptibility in Eq. ([5]), 



dct(l + go-zXnr) = 0. 



(30) 



The solutions of the real part of Eq. (|30[) give rise to the 
dispersions of the collective modes. Within the scheme of 
RPA, we find that the determinant det(l + ga z \nT) will 
develop a non-zero imaginary part when the dispersions 
fall in the particle-hole continuum, which is given in Fig. 
131 For this regime, the collective modes are unstable and 
will decay through particle-hole pairs. Substitute Eq. 
(|2"51) to Eq. ((50"]) . we obtain a gapless dispersion, which 
is given by the solution of 



(1 + 9Xnn)[l ~g(—+ y\nn)] + ^ffVxL = 0. (31) 



From Eq. (|28|) . we find that the density component is 
coupled with the transverse spin component. Therefore 
the zero sound is a coupled oscillation of the density and 
transverse spin excitations. There is a critical value for 
the interaction strength, i.e., g c = — , above which, the 
SOC Fermi gas experiences a Stoner instability. In the 
normal state regime we have g < g c . In Fig. 0] (a), we 
show the sound speed of this gapless mode by solving Eq. 
pip numerically. For comparison, the sound speed of the 
ordinary Fermi liquid without SOC is shown together. 

Substitute Eq. (|2T))) to Eq. (|3T)]) . we obtain one branch 
of gapped mode corresponding to transverse spin oscilla- 
tion. The frequency at |q| = of this mode is given by 
the solution of 



mg 



7T l&^ 2 e 2 F — oj* 
which yields the result 



= 0, 



At = 4\/l — mg/2TT'yEF, 



(32) 



(33) 



where is the energy gap of the transverse spin mode, 
mg = Eint/£kin is the dimensionless ratio of interaction 
energy and kinetic energy in 2D, which denotes the signif- 
icance of interaction strength in ultracold atomic system. 

Similarly, we consider susceptibility in the s L -s z sub- 
space in the limit q <^.kn and 7 < 1. The RPA suscep- 
tibility Xlz A 1S positive definite in the limit q — > while 
keeping y fixed in the normal state regime. Therefore 
there does not exist gapless excitations in the the lon- 
gitudinal and perpendicular spin subspace. In the limit 
q — > while keeping uj fixed, the interband susceptibility 
dominates and gives rise to 



m 

Xlz = — 



167 2 











(34) 



167 2 



The poles of RPA susceptibility Xlz A = XLz[^-gXLz] 1 
give rise to two branches of gapped modes corresponding 
to the longitudinal and perpendicular spin oscillations. 
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FIG. 5. (Color online) (a) The collective excitations for weak 
SOC with 7 = 0.01. The other parameters used here are the 
same with Fig. [4] The transverse, longitudinal and perpen- 
dicular spin excitations are labeled by T, L and Z. S denotes 
the zero sound mode, (b) The corresponding graphical solu- 
tions of the energy gaps at |q| =0 for 7 = 0.01. The energy 
gaps At and A_l for the gapped modes T and L are degen- 
erate, and the other solutions with nearly degenerated finite 
energy gaps labeled by G fall in the particle-hole continuum 
(the narrow pink region) and are therefore damped in weak 
SOC regime. 



FIG. 6. (Color online) The collective excitations for strong 
SOC with 7 = 0.5. The other parameters and the notations 
of these collective modes are the same with Fig. [5] We obtain 
one branch of gapless mode denoted as S and three branches 
of gapped modes denoted as T, L and Z. (b) The corre- 
sponding graphical solutions of the energy gaps at |q| = 
for strong SOC with 7 = 0.5. The energy gaps At and At 
are still degenerate with values close to the edge of the in- 
terband particle- hole continuum (shown in the inset figure). 
The other modes falling in the particle-hole continuum (the 
pink region) are damped. 



The frequencies at |q| = of these two modes are given 
as follows 

A L = 4-y/I - mg/2irje F7 (35a) 
A z = VI - mg/nje F , (35b) 

where Ar, and Az are the energy gaps for the longitudi- 
nal spin and perpendicular spin modes respectively. 

Compared with the formula of At in Eq. ([33]), we 
notice that the energy gaps for transverse spin and lon- 
gitudinal spin modes are degenerate: At — At- This 
degeneracy could be understood generally as follows. For 
finite momentum of the collective excitations q, the in- 
plane spin could be divided into transverse (s T ) and lon- 
gitudinal (s L ) components relative to the direction of q 
as defined in Sec. II. In the limit q — > 0, the in-plane 
spin components s T and s L are not distinguishable. And 
therefore the susceptibility for the in-plane spin fluctu- 
ations are isotropic. As a result the corresponding en- 
ergy gaps At and At equal to each other. For non-zero 
momentum, this degeneracy is lifted and the collective 
modes for s T and s L fluctuations have different disper- 
sions. This argument for the degeneracy also applies to 



the case for strong SOC, which is shown explicitly in the 
following discussions. 

In Fig. [SI we show the numerical results for collective 
modes and the graphical solutions of the energy gaps for 
weak SOC with 7 = 0.01, where the gapless zero sound 
mode is labeled as S and the transverse, longitudinal and 
perpendicular spin excitations are labeled by T, L and Z 
respectively. In the graphical solutions shown in Fig. [5] 
(b), we observed several additional solutions with nearly 
degenerated finite energy gaps labeled by G. However, 
all of gapped modes G fall in the particle-hole continuum 
of the interband excitations and are therefore damped. 

We now want to consider the collective modes for 
strong SOC. In the present experiments [IJ, [HI], the 
typical experimental values for the dimensionless ratio 
7 ranges from about 0.5 to 1. Thus we need to consider 
the case for 7 ~ 0(1) based on Eq. ([23 ]) and J3BJ 

in the following. The same as the case of weak SOC, 
the contribution from Eq. (|26|) is found to be negligible 
at the low-5 regime. We give the graphical solutions for 
these collective excitations in Fig. [5] and the sound speed 
of the gapless mode with strong SOC has been shown in 
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Fig. [1(b). 

The transverse, longitudinal and perpendicular spin 
excitation T, L, Z still exist as shown in Fig. [5] and 
the analytical formula for the energy gap of perpendic- 
ular spin modes At, A^, and Az given by (1551) and 
(l3"5"j) provide qualitative approximations with corrections 
about C(7 2 ). The two in-plane modes T, L exist close 
to the edge of the interband particle-hole continuum, 
which is shown by the inset of Fig. |6] (b). For typical 
experiment parameters, the dimensionless SOC strength 
7 ~ 0.5 — 1 and the dimensionless interaction strength 
mg = 2.5. The numerical analysis shows that the en- 
ergy gaps for the corresponding gapped modes in units 
of e F reads as A T = A L = 3.IO7 - 0.75j 2 + C(7 3 ), 
Az = I.8I7 + O.357 2 + C(7 3 ). Besides, there are several 
additional poles of determinant of the RPA susceptibil- 
ity with finite energy gaps. Similar to the case of weak 
SOC, all of these modes fall in the interband particle- hole 
continuum (shown in Fig. [3j and are therefore damped. 
The sound speed of the gapless mode as a function of 
the strength of SOC is shown in Fig. |4] (b), which de- 
creases as the strength of SOC is increasing. In partic- 
ular, the sound speed for strong SOC can also be deter- 
mined by the solutions of Eq. (|3Tj! • In fact, Eq. (|28"]) 
holds both for weak and strong SOC in the limit |q| — > 
with mui / nkp\(\\ kept fixed, which originates from the 
usual continuity equation and has been explicitly calcu- 
lated above. Similar forms could be found in graphene 
sheets [47|. 



V. COMPARISONS AND DISCUSSIONS 

We have noticed some previous works focused on col- 
lective modes in various solid and ultracold atomic sys- 
tems, such as the two dimensional electron gas (2DEG) 
in semi-conductors [48T4531 ]. the A phase and B phase in 
Helium-3 [HJ, [55|. Furthermore, the surface states of 
3D topological insulators are known as the helical liq- 
uid, which can be considered as the system with only the 
Rashba Hamiltonian. The collective modes in this helical 
liquid with Coulomb interactions or the repulsive short- 
ranged interactions have been explored in Ref. (56|. As 
a comparison and conclusion, we list these main results 
for the various systems in Table. |H 

For Coulomb interaction which is listed in the first col- 
umn of Table. |H the collective modes have the same 
dispersion relations lo cc y/q in 2D, since the poles of 
the RPA susceptibility just depend on the bare density- 
density susceptibility Xnn, which are solved by 

1 + Vc(<z)x«n = 0, (36) 

where Vc(q) is the Coulomb interaction vertex in 2D. 

In contrast to the Coulomb interaction, the s-wave in- 
teraction has two dramatically different features. First, 
the force range is short. Secondly, the s-wave interaction 
can be decomposed into the density and spin channels 



[30(. The interaction is repulsive in the density channel 
while attractive in the spin channel. Both of these two 
aspects have significant influences on the low-q behaviors 
of the collective modes. The results are listed in the sec- 
ond column of Table U which are classified as following: 

(1) . Ordinary Fermi liquid — In the absence of SOC, 
the bare density and spin susceptibility for the ordinary 
Fermi liquid is Xnnlix4- The fluctuation of spin chan- 
nel is suppressed because the interaction vertex for the 
spin susceptibility Xss is negative. The RPA equation 
1 + 9Xnn — gives rise to a linear dispersion relation, 
with the sound velocity v s > vp. The short range of the 
s-wave interaction is essential for the excitation of this 
mode. 

(2) . Helical liquid — Similarly, the collective modes in 
the helical liquid with short-ranged interaction are also 
gapless. However, in the helical liquid, the zero sound can 
be damped even for a repulsive interaction which is larger 
than the critical value U c as shown in Table Q] In fact, 
due to the a operator identity, j(x) = w/S(x) x z, which 
relates the density current to the in-plane component of 
the spin on the surface (56j , the spin and density fluctua- 
tions are intrinsically coupled. The effectively attractive 
interaction in the spin channel results in the damping of 
sound mode. 

(3) . SOC Fermi gas— The Hamiltonian of SOC Fermi 
gas includes both the quadratic term k 2 /2m and the 
Rashba SOC term. The Fermi momentum kp and recoil 
momentum fc# provide two scales for the system. The sig- 
nificance of the SOC is characterized by the dimensionless 
ratio 7 = kpi/kp. When we consider the gapped modes 
with energy gaps of 0{p/ep), the contribution of inter- 
band dominates while the intraband is well suppressed. 
We obtain three branches of undamped gapped modes 
bellow the interband particle-hole continuum correspond- 
ing to the spin oscillations (see Fig. [5] and |5J) . For the 
limit 7—^0, the gapped modes T, L, Z bellow the inter- 
band continuum will fall in the particle-hole continuum 
and are therefore damped. While for strong SOC limit 
7 — > 00, the energy gaps are much higher than the Fermi 
energy and are not observable in the low energy scales. 
These two limits coincide with the results of the ordi- 
nary Fermi liquid and the helical liquid. Besides, we also 
obtained one branch of gapless collective excitation, i.e., 
the zero sound with linear dispersion uo cx Vpq, which is 
a combined oscillation of density and transverse spin. 



VI. EXPERIMENTAL SIGNATURES AND 
SUMMARIES 

We have shown the behaviors of collective modes in the 
long wave length limit with SOC and repulsive s-wave 
interaction. Recently, the SOC in Fermi gas has been 
realized wtih 40 K atoms [l4[ and 6 Li atoms |3~5j] . In their 
experiments, the equal weight combination of Rashba- 
type and Dresselhaus-type SOC is realized. This is the 
first step towards the realization of pure type of SOC 
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TABLE I. Comparisons of the ordinary Fermi liquid, SOC Fermi liquid, and helical liquid. Both RPA results for the charged 
cases and the neutral cases with s-wave interaction are listed. 



Coulomb interaction 
=* 



s-wave interaction 



ordinary Fermi liquid lu 



/ 2nne 2 q 

me 



LU = ' VpQ 

^(l+Tr/mfl) 2 -! 

i): gapless w — v s q,v s > v F . 



SOC Fermi liquid 



Helical liquid 



2Vq 



ii): gapped 



A T = A L = 40 - mg/2Ti 1 e F , 



10 = n 



q 



lu = v F q{l + 8g 2 /g*),g/g c 0, 



* Here e is the dielectric constant of the ordinary Fermi liquid, and the e in the second and third row is the 
dielectric constant of semiconductors with SOC and the bulk of the topological insulator, respectively, n is the 
density of the charged particle |45j . 



Here V is the Drude weight, which is connected with the spin susceptibility [41 



experimentally. In this work, we focus on the Rashba- 
type SOC, and the Dresselhaus-type SOC is presented 
in Appendix B, which is demonstrated to give the same 
results as the Rashba SOC for collective behaviors. The 
short-ranged repulsive s-wave interaction can be achieved 
on the upper branch of a Feshbach resonance, where there 
are uncondensed Fermi gases in the absence of molecule 
formation Hi . The repulsive Fermi gas is metastable 
for observation when it is far away from the resonant 
regime, and has been successfully reached in the recent 
experiment (l7| . 

We choose the following typical experimental param- 
eters for quasi-2D system considered here. We consider 
about 10 4 40 K atoms in a pancake-shaped harmonic po- 
tential with the trapping frequencies 2tt x (10, 10, 400)Hz 
along the (x, y, z) direction. The system size is esti- 
mated as (37.8, 37.8, 5.98)/xm. The strength of SOC 7 
is chosen as 0.5 and a s is tuned to 2.25ao (within the 
normal state regime) with Feshbach resonance. For the 
parameters used here, we estimate that the zero sound 
velocity is about 1.09up, and the energy gaps for the 
gapped modes are At = Al = 1.23e_F and Az — 0.86e F . 
The Fermi velocity v F is about 0.028m/s and the Fermi 
energy e F — h x 33.36kHz. To observe these dynami- 
cal oscillations in experiment, we suggest the following 
methods, (i): The gapless mode could be excited by a 
short laser pulse focused near the center of the trap [57[ . 
The oscillation could be detected via the spatially re- 
solved images of the coupled density and transverse spin 
perturbances propagating through the trapped atomic 
cloud [531 . (ii): The gapped modes are spin oscillations, 
which are actually the oscillations of the internal hyper- 
fine states of atoms. The oscillations could be excited via 
a two-photon drive, and traced out through the state- 
selective absorption imaging method and repeating the 
experiment for many values of evolution time [HI, [59j] . 



In summary, we studied the SOC effect on the low-q 
behaviors of collective modes in Fermi gas with repulsive 
s-wave interaction. There are two categories of collective 
modes in this system. One branch has gapless dispersion, 
known as the zero sound. In SOC system, the density os- 
cillation is coupled with the spin oscillation for this mode. 
The other three branches are collective excitations with 
finite energy gaps. We calculated the sound speed of the 
gapless mode and the energy gaps of the gapped modes, 
and also estimated their values for typical experimental 
parameters. In addition, the repulsive s-wave interaction 
has a profound significance. In the presence of Coulomb 
interaction, the collective modes for ordinary Fermi liq- 
uid, SOC Fermi gas, and helical liquid have the same 
dispersion relation as lu ~ yfq in 2D. In contrast, the 
s-wave interaction leads to fundamentally different phe- 
nomena, such as the presence of linear dispersion of the 
zero sound and the gapped modes. The study on the 
collective modes of the SOC repulsive Fermi gas indi- 
cates some novel behaviors due to the presence of SOC, 
and also might have the immediate applicability to ex- 
perimental study of the SOC Fermi gases in the upper 
branch of the energy spectrum. 
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Appendix A: EVALUATIONS OF SOME 
RELEVANT INTEGRALS 



When we calculated the density and spin susceptibility 
in Sec. Ill, we met with some integrals of the azimuthal 
angle 9, such as 7o, 1%, I2, I3, 14, h- Since they have the 
similar structures, we evaluate them in this Appendix. 
The definitions of these integrations are listed as bellow: 



h(y) 
h(y) 
h(y) 
h(y) 
h(y) 



d9 



1 



2w y — cos 9 + i0 + ' 
dO cos 9 
2n y — cos 9 + i0 + ' 
dO cos 2 9 



2tt y ■ 
d9 



cos9 + i0+ , 
cos 3 9 



2tt y — cos ( 
d9 



i0+' 



sin 



2ir y — cos 9 + i0 + ' 



d6 



2-k y — cos 9 + i0 + 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



At first we note that if \y\ < 1, these integrals have a non- 
zero imaginary part; if \y\ > 1, these integrals are real. 
This can be seen from the following calculation. This 
property results in the damping of the collective modes 
in the particle-hole excitation continuum. 





(«):b[>l 



(b):\y\<l 



FIG. 7. Schematic of the integral path and the poles of the 
integrand functions. The location of the poles are different 
for \y\ > 1 and \y\ < 1, which are shown in (a) and (b) 
respectively. 

These integrals can be mapped into the integrals in 
the complex- z plane with a transformation z = e tS . The 
integral path C is the unit circle with the center at (0,0). 
We show the integral path and gives a schematic of the 
poles of the integrand function for \y\ > 1 and \y\ < 1 in 
Fig. [71 With the theorem of residue, these integrals are 



evaluated straightforwardly: 
sign(y) 



io(y) 



h{y) 



h{y) 



\fy rzr ^ 



e(M-i) 



? e(i-M), 



12/1 



iy 



e(M-i) 



e(i-M), 



y\y\ 



Vy 2 - 1 



o(|y|-i) 



=©(i -M), 



h{y) = -^-y 2 



w 



y 2 \y\ 
Vy^^i 



e(lwl-i) 



=0(1-12/1), 



h(y) =y- sign(y) - 19(\y\ - 1) 
h(y) 



- y 2 ^(i - \y\)}, 
4 + y 2 -lylV7 _ Te(|y|-i) 



- iyy/l ~ j/ 2 9 (1 - |y|). 



(A7) 

(A8) 

(A9) 

(A10) 
(All) 
(A12) 



where Q(x) is the unitstep function. When we were ex- 
ploring the analytical expressions of the susceptibility in 
Eq. (pD|). ([23} and (J26]) in the region q < k R , we used 
the asymptotic behaviors at \y\ — > 00: 



Jo (y)^-,Ii {y)^—2,I 2 (y) ~ — 



1 

2»' 



/,!,.^— ./,(</) ^^,J 5 ( y )c ^5. (A13) 



Appendix B: EQUIVALENCE BETWEEN THE 
RASHBA SOC AND DRESSELHAUS SOC 



We start with the single particle Hamiltonian wtih 
Dresselhaus SOC 111 



Hd = ^ + a {-ky<J x - k x a y ) - fi. (Bl) 

The various quantities calculated in this paper is closely 
based on the well-defined Feynman rules. The Feynman 
rules include the single particle Green's function and the 
interaction vertex. Within the scheme of RPA, the col- 
lective modes are determined by the poles of the RPA 
susceptibility, which is given by ([5]). It is apparent that 
the properties of the collective modes depends on the 
bare susceptibility x (k, w). Therefore, we want to derive 
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the relationships between the two types of SOC system 
as follows. 

For the Dresselhaus-type SOC, the non-interacting 
Green's function and the interaction vertex reads 

g d («, a,,) = E — /Vn+ - ( B2 ) 

w - £k. s + z0+ 

v£. rr ,(k,p,q) = ff/£(^A^)/,v(0p,0p-q)-(B3) 

where the upscript D represents the Dresselhaus-type 
SOC, and P® is the projection operator defined as 

P S D = i[l + s(- V« - fex^)]. (B4) 

The energy spectrum of the Dresselhaus-type SOC is the 
same with the Rashba-type. However, the spin polariza- 
tion is different from Rashba SOC, which is 

P* = -[1 + s(-k v a x + k x <i y )}. (B5) 

The trace in Eq. (flQ|) includes an overlap factor for 



the Dresselhaus case as 

F°<^ = tr [P S D (k - q/2) tjfjf (k + q/2) a v ] . (B6) 

Given Eq. (|B4j) and (|B5|) . the overlap factor for the two 
cases can be related by transformation k x — > —k Xl q x — > 
— (fc, while keeping the y component invariant. The bare 
susceptibility with Dresselhaus SOC is given by 

d.^lv i \ _ _ pg.M^ Z (^ k -q/ 2 - 6 ) ~ t (^ k +q/ 2 ^) 

(B7) 

Due to the rotation symmetry, we choose q = qe y for 
simplicity. Given the relationship between the overlap 
factor for Rashba SOC and Dresselhaus SOC and the 
rotational symmetry of the energy spectrum, we have 

X D ^(qe y ,uj) = X R ^( q e yi u J ). (B8) 

Therefore we conclude that all the properties and clas- 
sification of the collective modes for the two type of SOC 
are all the same exactly. 
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